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ABSTRACT 

In the subclass of high-mass X-ray binaries known as "microquasars" , relativistic hadrons in the jets 
launched by the compact object can interact with cold protons from the star's radiatively driven wind, 
producing pions that then quickly decay into gamma rays. Since the resulting gamma-ray emissivity 
depends on the target density, the detection of rapid variability in microquasars with GLAST and 
the new generation of Cherenkov imaging arrays could be used to probe the clumped structure of the 
stellar wind. We show here that the fluctuation in gamma rays can be modeled using a "porosity 
length" formalism, usually applied to characterize clumping effects. In particular, for a porosity length 
defined hy h = £/ f, i.e. as the ratio of the characteristic size £ of clumps to their volume filling factor 
/, we find that the relative fluctuation in gamma-ray emission in a binary with orbital separation a 
scales as ^Jhj-Ka in the "thin-jet" limit, and is reduced by a factor + (t>a/2£ for a jet with a 

finite opening angle <f>. For a thin jet and quite moderate porosity length h w 0.03 a, this implies a 
ca. 10% variation in the gamma-ray emission. Moreover, the illumination of individual large clumps 
might result in isolated fiares, as has been recently observed in some massive gamma-ray binaries. 
Subject headings: stars: binaries ~ stars: winds - gamma-rays: theory 



1. INTRODUCTION 

One of the most exciting achievements of high-energy 
astronomy in recent years has been to establish that high- 
mass X-ray binaries (HMXBs) and microquasars are vari- 
able gamma-ray sources (Aharonian et al. 2005, 2006; 
Albert et al. 2006, 2007). The variabihty is modulated 
with the orbital period, but in addition short-timescale 
flares seem to be present (Albert et al. 2007, Paredes 
2008). Since at least some of the massive gamma-ray bi- 
naries are known to have jets, interactions of relativis- 
tic particles with the stellar wind of the hot primary 
star seem unavoidable (Romero et al. 2003). At the 
same time, there are increasing reasons to think that 
the winds of hot stars have a clumped structure (e.g. 
Dessart & Owocki 2003, 2005; Puis et al. 2006, and ref- 
erences therein). The observational signatures of such 
clumping often just depend on the overall volume flU- 
ing factor, with not much sensitivity to their scale. Here 
we argue that gamma-ray astronomy can provide new 
constraints on the clumped structure of stellar winds in 
massive binaries with jets. At the same time, our anal- 
ysis provides a simple formalism for understanding the 
rapid flares and flickering in the light curves of these ob- 
jects. Our basic hypothesis is that the jet produced close 
to the compact object in a microquasar will interact with 
the stellar wind, producing gamma-rays through inelas- 
tic pp interactions, and that the emerging gamma-ray 
emission will present a variability that is related to the 
structure of the wind. Thus the detection of rapid vari- 
ability by satellites like GLAST and by Cherenkov arrays 
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like MAGIC II, HESS II, and VERITAS can be used as 
a diagnosis of the structure of the wind itself. 

2. JET-CLUMP INTERACTIONS 
2.1. The general scenario 

The basic scenario explored in this paper is illustrated 
in figure[TJ A binary system consists of a compact object 
(e.g., a black hole) and a massive, hot star. The com- 
pact object accretes from the star and produces two jets. 
For simplicity, we assume that these jets are normal to 
the orbital plane and the accretion disk (see otherwise 
Romero & Orellana 2005). We also assume a circular 
orbit of radius a. The wind of the star has a clumped 
structure and individual clumps interact with the jet at 
different altitudes, forming an angle ^' with the orbital 
plane. The z-axis is taken along the jet, forming an angle 
9 with the line of sight, with the orbit in the a;j/-plane. 
The jet has an opening angle 0. To consider the effects 
of a single jet-clump interaction, we first adopt a model 
for the jet^ (Sect. W2\ . 

In addition to wind clumping, there can also be in- 
trinsic variability associated with the jet. This includes 
orbital modulation, as observed in LS 5039 or LS I 61 303 
(e.g., Aharonian et al. 2006; Albert et al. 2006), and pe- 
riodic precession of steady jets (e.g., Kaufman-Bernado 
et al. 2002). Both these long-term, periodic varia- 
tions would be quite distinct from the rapid, stochas- 
tic variations from wind clumps. Intrinsic disturbances 
and shocks in jets can produce aperiodic variability that 
might be confused with variability associated with jet- 
clump interactions. In microquasars such intrinsic fluc- 
tuations are expected to arise in the context of the jet- 
disk coupling hypothesis, as proposed by Falcke & Bier- 
mann (1995) for the case of AGNs, and observationally 

^ Note that the interaction of a beam of protons and a cloud 
from a star has been considered before, in the context of pulsars, 
by Aharonian & Atoyan (1996). An early report of some material 
presented here can be found in Romero et al. (2007). 
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Fig. 1. — Sketch of the assumed model, described further in the 
text. 

demonstrated for a galactic microquasar by Mirabel et 
al. (1998). The same effect has been observed in AGNs 
(Marscher et al. 2002). Thus intrinsic variability in the 
jet would likely be preceded by a change in the accre- 
tion disk X-ray activity, whereas in the case of a jet- 
clump interaction, the effect should be the opposite: first 
the gamma-ray flare would appear, and then, an non- 
thermal X-ray flare produced by the secondary electrons 
and positrons as well as the primary electrons injected 
into the clump would show up. Depending on the mag- 
netic field and the clump density, the X-ray radiation 
could be dominated by synchrotron, inverse-Compton, or 
Bremsstrahlung emission, with a total luminosity related 
to that of the gamma-ray flare. In summary, simulta- 
neous X-ray observations with gamma-ray observations 
could be used to differenciate jet-clump events from in- 
trinsic variability produced by the propagation of shocks 
in the jets. 

2.2. Basics of the jet model and jet- clump interaction 

The matter content of the jets produced by micro- 
quasars is not well-known. However, the presence of 
relativistic hadrons in the jets of SS433 has been di- 
rectly inferred from iron X-ray line observations (e.g. 
Kotani et al. 1994, 1996; Mighari et al. 2002). In ad- 
dition, the large perturbations some jets cause in the in- 
terstellar medium imply a significant baryon load (Gallo 
et al. 2005, Heinz 2006). The fact that the jets are 
usually well-coUimated also favors a content with cold 
protons that provide confinement to the relativistic gas. 
We adopt here the basic jet model proposed by Bosch- 
Ramon et al. (2006), where the jet is dynamically dom- 
inated by cold protons. Since the jet launching likely 
stems from magneto-centrifugal effects (e.g., Blandford 
& Payne 1982), the jet magnetic field is assumed to be 
in equipartition with the particle energy density, with 
typical values of 1 kG. 

Shocks from plasma collisions in the jet can produce a 
non-thermal relativistic particle population. But only 



a fraction qj w 0.1 of the total jet luminosity Lj w 
lO'^^ergs"^ is expected to be converted into relativistic 
protons by such diffusive shock acceleration at the jet 
base (e.g., Riegler ct al. 2007). The resulting gamma-ray 
emission can be calculated as in Romero et al. (2003) and 
Orellana et al. (2007). For interaction between relativis- 
tic (~TeV) protons in the jet with cold protons in the 
wind, a characteristic cross section iso'R:i3.4xlO~^^ cm^ 
(Kelner et al. 2006). For a typical wind mass loss rate 
M = lO"^M0yr-i and speed v = 1000kms"\ the char- 
acteristic wind column depth traversed by the jet from 
an orbital separation distance a ~ 0.2AU ~ 3 x 10^^ cm 
is TV « 5 X 10^^ cm~^. This implies that only a small 
fraction, ~ crN « 0.002, of relativistic particles in the 
jet are converted to gamma-rays by interaction with the 
entire wind. The leads to a mean gamma-ray luminosity 
of — qjT-ujLj w 2 X lO^^ergs"^. 

Clumps in the wind can lead to variations and flares in 
this gamma-ray emission. For clumps of size £ ~ 10^^ cm, 
corresponding to a few percent of the stellar radius, 
the flow into the jet at the wind speed implies a flare 
timescale less than an hour. While quite short, this is 
comparable to the variability already detected in Cygnus 
X-1 by MAGIC (Albert et al. 2007). HESS II and 
MAGIC II will have a higher sensitivity, so these instru- 
ments should be able to detect variability from galactic 
sources like LS 5039, LS I +61 303 and Cygnus X-1 on 
timescales below an hour. 

The flare brightness depends on the clump column 
depth and the resulting fraction of the relativistic particle 
luminosity converted to gamma rays. For clumps of the 
above size with a volume filling factor / = 0.1, using the 
above wind parameters at an orbital separation distance 
gives a clump column iVc = 3 x lO^^cm^^. The asso- 
ciated clump attenuation fraction is = crNc ~ 10^**, 
implying a flare gamma-ray brightness of L-^ — qjTcLj w 
lO'^^ergs"^. Stronger flares could result when a large 
clump crosses close to the base of the jet. Overall, if a 
microquasar is observed in an active state (i.e. when the 
jet is powerful), then satellite instruments like GLAST 
and ground-based Cherenkov telescopes should be able 
to detect variability down to timescales of ^ 1 h or so, 
sufficient to measure variations associated with jet inter- 
actions with wind clumps. 

The above picture assumes that the jet is not signifi- 
cantly dispersed or attenuated by other interactions with 
the wind, for example gyro-scattering off magnetic field 
fluctuations in the clumps. Taking a characteristic tem- 
perature T « 10* K along with the above parameters for 
the wind and clumps, we can estimate that at the orbital 
separation distance a = 0.2 AU, clumps have a typical 
thermal energy density i?c = (3/2)nc/cT ~ 0.07 erg cm""^, 
with the corresponding equipartition magnetic field thus 
of order a Gauss. For relativistic protons of Lorcntz fac- 
tor 7, the associated gyroradius is just Vg = 3O7 km. 
Even for TeV particles with 7 « 10^, this is much less 
than the clump size, rg <^ £, implying that individual 
such particles should be quite effectively deflected by 
such clumps. 

However, this does not mean that such gyroscattering 
by wind clumps can substantially disperse the jet. The 
simple reason is that the energy density of the jet com- 
pletely overwhelms that of the wind clumps. For a jet 
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with opening 



1° and thus soHd angle fl = iKfP' 



10"-^ stcr, the energy density at an orbital distance a 
is Ej = Lj/(p.ca^) « 4 X 10'*ergcm~^, nearly a million 
times higher than for the clumps. This suggests that, 
while clump-jet interactions may substantially perturb 
or even destroy the clumps, the back-effect on the jet 
should be very small. Moreover, while the dynamics of 
such clump destructions are likely to be complex, the 
overall exposure of clumped wind protons to interaction 
with the relativistic protons in the jet may, to a first 
approximation, remain relatively unaffected. Overall, it 
thus seems reasonable to assume a simple interception 
model of jet-wind interaction, withrelatively little dis- 
persal or attenuation of the jet through the wind. 

3. POROSITY-LENGTH SCALING OF GAMMA-RAY 
FLUCTUATION FROM MULTIPLE CLUMPS 

Individual jet-clump interactions should be observable 
only as rare, flaring events. But if the whole stellar 
wind is clumped, then integrated along the jet there will 
be clump interactions occurring all the time, leading to 
a flickering in the light curve, with the relative ampli- 
tude depending on the clump characteristics. Under the 
above scenario that the overall jet attenuation is small, 
both cumulatively and by individual clumps, the mean 
gamma-ray emission should depend on the mean num- 
ber of clumps intersected, while the relative fluctuation 
should (following standard statistics) scale with the in- 
verse square-root of this mean number. But, as we now 
demonstrate, this mean number itself scales with the 
same porosity-length parameter that has been used, for 
example, by Owocki and Cohen (2006) to characterize 
the effect of wind clumps on absorption of X-ray line 
emission (see also Oskinova, Hamman, and Feldmeier 
2006). 

Let us again consider the gamma-ray emission inte- 
grated along the jet. Representing the relativistic par- 
ticle component of the jet as a narrow beam with con- 
stant luminosity Lf, = qjLj along its length coordinate 0, 
the total mean gamma-ray luminosity scales (in the 
small-attenuation limit L-y Lb) as 



r 

{L^) = L^a / n'{z) dz , 
Jo 



(1) 



where n{z) is the local mean wind density (i.e. averaged 
over any small-scale clumped structure), and a is the 
gamma-ray conversion cross-section defined above. 

The fluctuation about this mean emission depends on 
the properties of any wind clumps. A simple model as- 
sumes a wind consisting entirely of clumps of charac- 
teristic length £ and volume filling factor /, for which 
the mean-free-path for any ray through the clumps is 
given by the porosity length h = £/f. For a local in- 
terval along the jet Az, the mean number of clumps in- 
tersected is thus ANc = Az/h, whereas the associated 
mean gamma-ray production is given by 



(ALj) — L^anAz = LbanANch. 



(2) 



But by standard statistics for finite contributions from a 
discrete number AN^, the variance of this emission about 

the mean is 

(AX^) - {AL,f = ^^^^^ = Lla-n-hAz . (3) 



Each clump-jet interaction is an independent process; 
thus, the variance of an ensemble of interactions is just 
the sum of the variances of the individual interactions. 
The total variance is then just the integral that results 
from summing these individual variances as one allows 
Az dz. Taking the square-root of this yields an ex- 
pression for the relative rms fluctuation of intensity. 



v?hdz 



(4) 



Note that, in this linearized analysis based on the weak- 
attenuation model for the jet, the cross-section a scales 
out of this fluctuation relative to the mean. 

As a simple example, for a wind with a constant veloc- 
ity and constant porosity length /i, the relative variation 
is just 




' dxl{\ 



dxjil^x^) 



ira . 



(5) 



Typically, if, say h « 0.03a, then SL^/L^ « 0.1. This im- 
plies an expected flickering at the level of 10% for a wind 
with such porosity parameters, occuring on a timescale 
of an hour or less. 

4. GAMMA-RAY FLUCTUATIONS FROM A FINITE-CONE 

JET 

Let us now generalize this analysis to take account of a 
small but finite opening angle 4> for the jet cone. The key 
is to consider now the total number of clumps intersecting 
the jet of solid angle f2 « At a given distance z from 
the black hole origin, the cone area is fiz^ = {(pz)^ . For 
clumps of size £ and mean separation L, the number of 
clumps intercepted by the volume Q.z'^Az is 



AiV^ = Az 



Az 



[1 + {(bz/£f 



(6) 



where the latter equality uses the definition of the poros- 
ity length h = £/ f in terms of clump size £ and volume 
filling factor / = £^/L^. 

Note that the term "intercepted" is chosen purpose- 
fully here, to be distinct from, e.g., "contained". As 
the jet area becomes small compared to the clump size, 
the average number of clumps contained in the vol- 
ume would fractionally approach zero, whereas the num- 
ber of clumps intercepted approaches the finite, thin-jet 
value, set by the number of porosity lengths h crossed 
in the thickness Az. As such, for (j)z <C £, this more- 
general expression naturally recovers the thin-jet scaling, 
AA^c = Az/h, used in the previous subsection. 

Applying now this more-general scaling, the emission 
variance of this layer is given by 



AA^, 



Lla'^n^hAz 



(7) 



Obtaining the total variance again by letting the sum be- 
come an integral, the relative rms fluctuation of intensity 
thus now has the corrected general form. 



5L, _^j!^n^hdz/{l + <l>H^/£^) 
{L^) n dz 



(8) 



4 



Owocki et al. 



For the simple example that both the porosity length 
h and clump size £ are fixed constants, the integral forms 
for the relative variation becomes 




dx/[{l+p^x^){l + X- 



2\2] 



(9) 



where p = (f)a/£ defines a "jet-to-clump" size parameter, 
evaluated at the binary separation radius a. Carrying 
out the integrals, we find the fluctuation from the thin- 
jet limit given above must now be corrected by a factor 



Cp — 



1 



1 



(10) 



where the latter simplification is accurate to within 6% 
over the full range of p. 

In the thin-jet limit p — (pa/f. <^ 1, the correction 
approaches unity, as required. But in the thick-]ei limit, 
it scales as 




l/a. 



(11) 



When combined with the above thin-jet results, the 
general scaling of the fluctuation takes the approximate 
overall form 

5Lj_^ j h/TTa 

{L^) ~ Y 1 + q^a/2£ ^ ' 

wherein the numerator represents the thin-jet scaling, 
while the denominator corrects for the finite jet size. 

If the jet has an opening of one degree, then (j) = 
(tt/ISO) « 1.7 X 10^^ radian. If we assume a clump 
filling factor of say, / = 1/10, then the example of the 
previous section for a fixed porosity length h = 0.03 a 
implies a clump size £ = 0.003 a, and so a moderately 
large jet-to-clump size ratio of p ~ 6. But even this gives 
only a quite modest reduction factor Cp ~ 0.5, yielding 
now a relative gamma-ray fluctuation of about 5%. 

The bottom line here is thus that the correction for 
finite cone size seems likely to give only a modest (typi- 
cally a factor two) reduction in the previously predicted 
gamma-ray fluctuation levels of order 10%. This holds 
for clump scales of order a few thousandths of the binary 
separation, and for jet cone angles of about 1 degree. As 
the ratio between these two parameters decreases (still 
keeping a fixed porosity length), the fluctuation level 
should decrease in proportion to the square root of that 
ratio, i.e. SL^ oc \JTf^ oc 1/ ^Jp. 



5. CONCLUSION 



Overall, for a given binary separation scale a, our gen- 
eral model for gamma-ray fluctuation due to jet interac- 
tion with clumped wind has just two free parameters, 
namely the porosity length ratio h/a, and the jet-to- 
clump size ratio p — <j)a/£. Given these parameters, 
then, within factors of order unity, the predicted rela- 
tive gamma-ray fluctuation is given by eqn. (|12p . For 
reasonable clump properties with h « 10^ « 0.03a, the 
fluctuation amplitude would be a few percent. 

Note however that the formalism here is based on a 
simple model in which all the wind mass is assumed to 
be contained in clumps of a single, common scale £, with 
the regions between the clumps effectively taken to be 
completely empty. More realistically, the wind structure 
can be expected to contain clumps with a range of length 
scales, superposed perhaps on the background smooth 
medium that contains some nonzero fraction of the wind 
mass. For such a medium, the level of gamma-ray fluc- 
tuation would likely be modified from that derived here, 
perhaps generally to a lower net level, but further anal- 
ysis and modeling will be required to quantify this. 

One potential approach might be to adopt the "power- 
law porosity" formalism developed to model the effect of 
such a clump distribution on continuum driven mass loss 
(Owocki, Gayley, and Shaviv 2004). This would intro- 
duce an additional dependence on the distribution power 
index ap, with smaller values ctp — > tending to the 
smooth flow limit. But for moderate power indices in 
the range 0.5 < ap < 1, we can anticipate that the above 
scalings should still roughly apply, with some reduction 
that depends on the power index a^, if one identifles the 
assumed porosity length h with the strongest clumps. 

Thus while there remains much further work to deter- 
mine the likely nature of wind clumping from hydrody- 
namical models, the basic porosity formalism developed 
here does seem a promising way to characterize its broad 
effect on key observational diagnostics, including the rel- 
ative level of fluctuation in the gamma-ray emission of 
HMXB microquasar systems. 
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